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1.  Introduction 

The  subject  of  how  a  boundary  layer  responds  to  fluctuations  in  the 
external  flow  about  a  steady  mean  has  been  discussed  recently  by  several 
authors  [6,  7,  8,  9,  10] .  One  reason  for  an  interest  in  such  problems  arises 
in  a  study  of  the  "Rijke  tube"  phenomenon  [  6,  7,  11] .  There  it  appears  to  be 
important  to  know  the  phase  lag  in  the  heat  transfer  from  a  heated  gauze 
situated  in  a  stream  with  an  unsteady  velocity.  Lighthill  [  7]  investigated 
the  oscillating  flow  problem  using  boundary  layer  theory  and  the  Karman- 
Pohlhausen  method.  It  is  apparent  that  his  results  can  not  be  valid  near  the 
leading  edge  of  the  object  (say  a  heated  ribbon).  However  most  of  the  heat 
release  takes  place  at  the  leading  edge;  hence  it  appears  desirable  to  study 
the  behaviour  of  the  velocity  fluctuations  near  the  leading  edge. 

In  this  paper  the  viscous  incompressible  flow  past  a  semi-infinite 
flat  plate  is  investigated  when  the  flow  at  infinity  is  UQ  +  exp  (iuit).  The 
analysis  is  based  on  a  system  of  equations  derived  by  a  modified  Oseen 

linearization  of  the  equations  of  motion  [1,2].  It  is  found  that  when 

y 

(wV)'1'  »  U0  and  Uj  the  flow  field  consists  of  the  sum  of  a  mean  flow 
(due  to  U0)  and  an  oscillatory  motion  (due  to  exp  (icot) | .  The  phase  advance 
in  the  oscillatory  motion  is  w/8  at  the  leading  edge  and  increases  monotonically 
to  ji/4  in  a  distance  approximately  given  by  (u>/2 v^)  lx  =  2.5.  The  result  for 

j/ 

(o;/2V)a  x  ^  2.5  is  a  shear  wave  type  solution  and  corresponds  to  an  infinite 
lamina  performing  harmonic  oscillations  in  a  viscous  fluid. 

In  Sec.  5  the  problem  is  considered  when  the  interaction  of  the  mean 
motion  and  the  oscillatory  motion  can  not  be  neglected.  Expressions  for  the 
amplitude  and  phase  advance  of  the  unsteady  skin  friction  near  the  leading 
edge  are  found.  In  particular  the  phase  advance  at  the  leading  edge  increases 
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monotonically  from  0  to  n/8  as  the  dimensionless  parameter  (wV)  Y(c  UQ) 
increases  from  0  to  oo.  Here  c  is  a  constant  which  is  defined  in  Sec.  2. 
Far  down  the  plate  the  solution  approaches  the  shear  wave  solution,, 

2.  Formulation  of  Problem 

The  equations  governing  a  two-dimensional  viscous,  incompressible 
flow  are 


ux  +  vy  =  0  » 

/(ut  +  uux  +  vuy)  =  -  px  +  /^Au, 

/(vt  +  UVX  -  VVy)  =  -  Py  +  ytaAV  . 


(2.1) 


Here  x  and  y  are  the  usual  Cartesian  coordinates,  u  and  v  are  the  corresponding 
velocities,  p  is  the  pressure,  J° t he  density  (which  is  taken  constant),  V  =  /*/f 
is  the  kinematic  viscosity,  A  is  the  Laplacian  operator,  and  x,  y,  or  t  used 
as  a  subscript  denotes  partial  differentiation  with  respect  to  that  variable. 

We  wish  to  determine  the  velocity  field  associated  with  a  semi -infinite  flat 
plate  situated  in  an  unsteady  stream.  For  convenience  we  shall  assume  that 
the  semi-infinite  flat  plate  occupies  the  half  line  y  =  0,  x  >,  0.  In  particular 
we  wish  to  consider  this  problem  when  the  imposed  velocity,  U(t),  far  from 
the  plate  is  of  the  form 

U(t)  =  u0  *  U,  e1"1  ,  (2.2) 


where  UQ  and  Uj  are  constants.  The  boundary  conditions  which  u  and  v  must 
satisfy  are 


u  =  0,  v  =  0  on  y  =  0,  x  ?  0, 
u— >U0  +  Ul  e^1,  v  — >  0 


(2.3) 


far  from  the  plate. 
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In  order  to  treat  this  problem  we  linearize  Eqs.  (2.1)  by  replacing 

the  convective  terms  uux*+  vuy ,  and  +  vvy  by  c  U  ux  and  c  U  vx.  Here  c 

is  a  constant  which  in  general  should  depend  upon  the  parameters  involved 

in  the  problem.  The  case  Uj  =  0,  c  =  1  gives  the  well  known  Oseen 

approximation.  The  modification  of  the  Oseen  method  (i.e.  c  f-  1)  has  been 

used  in  treating  the  steady  flow  over  a  semi -infinite  plate  [  1] .  Also  it  has 

problems 

been  applied  to  several  other  steady  flov^,  and  its  rationalization  discussed 
in  [  2] .  In  [  2]  it  was  noted  that  for  small  Reynolds  numbers  that  the  value 
c  =  .43  gave  results  in  good  agreement  with  experiment  in  all  cases  considered. 
When  the  flow  is  unsteady,  the  choice  of  c  is  clearly  more  complicated. 

However  in  a  recent  publication  [  3]  this  type  of  linearization  was  used  to 
compute  the  fluctuating  heat  release  from  a  heated  ribbon  when  U(t)  had  the 
form  given  by  Eq.  (2.2),  and  the  choice  of  c  =  .43  gave  excellent  results. 

The  form  of  U(t)  suggests  that  we  write  u  and  v  as 

u(x,  y,  t.)  =  uQ  (x,  y)  +  ut  (x,  y)  e^1  , 

=  U0|lTXy(x,  y)}+U1{l  +  ^y(x,y)}elu)t  ) 
and  ^(2.4) 

V(X,  y,  t)  =  vQ(x,  y)  +  v,  (x,  y)  emt  , 

=  -  U  X  (x,  y)  -  U,  ^  (x,  y)  e1"* . 

U  X  A  . _ / 

In  order  to  obtain  a  solution  of  this  form  it  is  necessary  to  neglect  quadratic 
terms  in  ux ;  i.e.  we  must  restrict  ourselves  to  cases  in  which  Uj  is  small 

y 

compared  to  either  UQ  or  (cov1)4'.  Then  using  the  linearization  discussed 
above  and  Eqs.  (2.4)  we  obtain  the  following  equations  for  X  and 

AA  X-  aAX  =  0  ,  (2.5) 

A 


AAf  -  aAfx  -  /32  A^=  a^Xx  . 


(2.6) 
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Here  a  =  cUq/v  ,  and  fi*  =  ico/vJ  .  The  boundary  conditions  that  Pi  and  [f 
must  satisfy  are 


V°<  ^x=°;  "  *  *>  V  -1 

o,  Xy-^o,  yx-+  0,  yy— < 


on  y  =  0,  0<  x*co, 


0, 


(2  c  7) 


far  from  the  plate.  The  conditions  and  Xx  =  0  on  y  =  0,  0*  x  *  oo  clearly 
hold  for  all  x  and  we  satisfy  them  by  requiring  V''  and  X  be  constant  on  y  =  0: 
in  particular  Y'  (x,  o)  =  X  (x,  o)  =  0,  As  we  shall  see,  when  u>  is  very  large 
the  terms  involving  a  in  Eq„  (2,6)  may  be  neglected  and  we  obtain  homogeneous, 
independent  equations  for  X  and  Y/.  That  is  to  say,  the  velocity  distribution  is 

simply  the  flow  due  to  a  velocity  U0  at  infinity  plus  that  due  to  a  velocity  eiu' 

at  infinity.  In  general  of  course  these  terms  cannot  be  neglected,  and  hence  the 
interaction  of  the  mean  flow  and  the  time  dependent  flow  must  be  considered, 

The  term  aAyx  represents  the  convective  effect  of  the  mean  flow  on  the  time 
dependent  flow,  while  the  non -homogeneous  term  aA")^  represents  the  convective 
effect  of  the  time  dependent  flow  on  the  mean  flow. 

To  complete  the  formulation  of  our  problem  we  must  note  that  while 
X,  ■/,  Xy,  Yy,  Xyyy  and  Yyyy  are  continuous  across  the  plate  Xyy  and  yr 
which  are  proportional  to  the  velocity  gradient  are  discontinuous.  We  may 
express  this  discontinuity  as 

Xyy(x,  0+)  =  “  *yy(x’  °’)  =  2  f°(x)  ’ 

Y'yy (x,  0+)  =  *  tyy(x,  0  *)  =  |  fi  (x)  . 


The  functions  fQ(x)  and  ft  (x)  are  identically  zero  for  0.  Now  the  skin 
friction,  V  ,  is  given  by 


r  =/4 


<)u 


y=o 


v. 


~W 


,iu>t 


y=  0 
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Upon  substituting  for  u0  and  ux  from  Eqs,  (2.4)  and  using  the  above  relations 
we  may  express  V  m  the  form 


r  *  r0  +  r,  e 


icot 


V*uo 


Ux)  aiuot 


(2.8) 


In  what  is  to  follow  we  shall  be  primarily  concerned  with  the  determination 
of  fQ(x)  and  fj  (x)„ 

3.  Case  U,  =  0 

In  this  case  our  problem  reduces  to  a  solution  of  the  problem  of  the 
flow  over  a  semi-infinite  flat  plate  due  to  a  flow  U0  at  infinity.  Thus  we  must 
solve  Eq.  (2.5)  subject  to  the  boundary  conditions  (2.7).  Even  when  f  0  we 
still  need  the  solution  of  this  problem  in  order  to  obtain  ^  (see  Eq.  (2.6) ). 

The  solution  of  this  problem  has  been  given  in  [  1] .  Since  the  analysis  of  the 
general  problem  will  follow  the  lines  of  [  1] ,  but  is  much  more  complicated 
in  its  detail;  and  also  since  we  need  the  results  found  in  [  1]  it  is  convenient 
to  use  this  simpler  problem  to  illustrate  the  method  of  analysis  that  is  to  be 
used.  Hence  we  shall  repeat  briefly  part  of  the  computations  described  in  [  1] . 

Our  problem  can  be  most  conveniently  treated  by  the  use  of  Fourier 
transforms.  We  define 

••  O*  Qg 

Xa,1!)-//  e'i(}x+lyh(x, y)dxdy,  TQ(f )  =  /Vifx  fQ(x)  dx. 

(3.1) 

In  all  that  is  to  follow  a  barred  quantity  will  denote  the  Fourier  transform 
of  that  quantity.  Taking  the  Fourier  transform  of  Eq.  (2.5),  and  using  the 
boundary  conditions  and  discontinuity  relation  for  X  we  obtain 


(?2  +  ?2)<5!+f  Tla?)X(f,?)  =  i710(f)  . 


(3.2) 
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In  order  to  avoid  difficulties  later  it  is  desirable  to  write  this  equation  as 
the  limit  as  k — ^  0  of  the  equation4 


+  f  +  k2}  +(f,  ia)  (J-  ik)}X  (5,?)  =  i?£0fi).  (3.3) 


If  we  let  X(^,  yf  denote  the  inverse  transform  of  X  (?,^)  over  ^  we  obtain 
upon  solving  Eq.  (3.3)  and  using  simple  contour  integration  that 


ei^yX(F,'rl)  d^  =  J  e_i^xX(x,y)  dx  , 


(3.4a) 


exp{-jy|  (f  +  k2)*}-exp{-|y/(f>  ia.f  (l~  ik)*}  _ 


2i(a  -  k)  (?-  ik) 


*n(?) 


o 

(3,4b) 


In  particular  we  obtain  from  Eq.  (3.4b)  that 


Xy(F,  y  =  0)  = 


-  fo(^)/2 

(F*  ik)*  ■  (?+  ia)*-»-  (£+  ik)*} 


(3  5) 


A  second  equation  for  Xy(£,  y  =  0)  can  be  obtained  by  using  our 
boundary  conditions.  First  however  the  condition  that  Xy  =  -  1  on  the  plate 
must  be  written  as  the  limit,  as  a — >  0,  of  Xy(x,  o)  =  -  exp  (-ax)  on  x  ?  0. 
(The  remarks  made  about  k  previously  also  apply  to  of.)  Now  differentiating 
(3.4a)  with  respect  to  y  and  setting  y  =  0  we  have 


*y  (f,  y  =  o)  =  j  el5x/y(x,  y  =  0)  dx  , 


=  V(?) 


1 

!(?-  1«)  ’ 


(3.6) 


where 


v(?)  = 


Xy(x,  y  =  0)  dx  . 


+Here  and  throughout  the  rest  of  this  paper  the  symbol  k  will  be  used  to 
represent  a  small,  positive,  real  quantity  which  is  to  be  shrunk  to  zero 
at  an  appropriate  moment. 

*We  shall  be  careful  to  write  X(£,  y)  if  we  mean  the  x  transform  of  X(x,  y) 
and  X(S\  *?)  if  we  are  referring  to  the  two-dimensional  transform. 
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The  function  v(  f)  is  of  course  not  known.  Equating  (3.6)  and  (3.5)  we  obtain 


vU)  -  - i - 

i(?  -  ia) 


-fq(S)/2 

( ?-  ik)1  ik)'£  .  (5+  ia)* 


(3.7) 


We  can  now  determine  v(  ?)  and  fQ(  ?)/2  by  the  Wiener -Hopf  technique. 

If  we  anticipate  that  Xy(x,  y  =  0)  decays  exponentially  in  (-x)  on  the  negative 
real  axis,  then  v7 (?)  is  an  analytic  function  of  the  complex  variable  ?  in  some 
upper  half -plane  (UHP)  which  includes  the  real  axi£.  Similarly  f0(?)  is  analytic 
in  some  lower  half -plane  (LHP)  including  the  real  axis.  Also  l/(?  -  ia)  is 
analytic  in  the  LHP,  Im(?)<a.  And  finally  K(  ?)  =  (f  -  ikf  |(f  +  ik)x  +  (?+  ia)*j“ 
is  analytic  in  the  strip  -  k  <  Im(? )  <  k.  Let  us  assume  that  we  can  split  Kf?) 
into  the  product  of  two  factors  K  (1)  and  K.(?)  such  that  K+(?)  is  analytic  in  the 
UHP,  Im(?)  >  -  k;  and  K.C?)  is  analytic  in  the  LHP  Im(]j)^  k.  Then  we  may 
rewrite  Eq.  (3.7)  as 


f0(?)/2  K  (ia) 

^  i  i  n  —  sm  ....  i  ,TTi— 

K.(?)  i(?  -  ia) 


*  v(  ?)  K+(f )  - 


K+(ia)  -K+(?) 
i(f  -  ia) 


=  E(?) 


(3.8) 


The  left-hand  of  this  equation  is  analytic  in  the  LHP,  while  the  right  side  is 
analytic  in  an  overlapping  UHP,  and  hence  they  define  an  entire  function  EC?). 
By  using  order  conditions  of  the  known  functions  at  infinity  and  the  fact  that 
E(?)  is  an  entire  function  we  can  show  that  E(f )  =0,  hence 

Ml)  .  K-<?>  (3.9) 

2  i(|  -  lor) 

We  recall  once  again  that  we  are  interested  in  evaluating  f0(^)/2  in  the  limit 
a  — >  0,  k  — >  0. 

In  this  particular  problem  K_(  ?)  and  K+(0)  can  be  determined  by 
inspection  of  K(f)  (we  will  not  be  so  fortunate  later  on).  We  see  that 
K.(?)  =  (X  -  ik)*  and  K+(0)  =  (ia)*  in  the  limit  as  a  and  k — >0.  Substituting 
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in  (3.9)  we  may  write  fQ(f)  as 


U!)/2  =  ->TT7*~ 
0  (it  +  k)* 


(3.10) 


We  have  written  (3 <,10)  in  the  preceding  manner  to  indicate  clearly  that  the 
singularity  is  in  the  UHP„  In  the  limit  as  k — we  find  by  the  conventional 
inversion  formula  that 


yx)/2  = 


a  ^ 
nx 


0 


cU0Y* 


7t  V>  X 


x  0 


x  ^  0 


(3.11) 


As  mentioned  previously  the  choice  of  c  =  .43  gives  good  agreement  with  the 
non-linear  theory  and  with  experimental  results  for  small  Reynolds  number. 
However  in  this  case  the  Reynolds  number  is  infinite  and  the  result  which  is 
consistent  with  the  Blasius  solution  is  given  by  c  =  .35  [  1] .  It  is  clear  that 
the  numerical  choice  of  c  can  not  affect  the  character  of  the  result,  and  hence 
for  any  value  of  c  in  this  range  the  result  will  describe  the  basic  features  of 
the  flow.  Finally  knowing  f0(?)  it  is  possible  to  invert  X(?  ,  y)  and  obtain 


X<s,y)=^  |2(r+x)* 


(2  jra)^ 


1  -  e 


■a(x  -r)/2 


y  erfc 


a(r  -  x) 


where  r2  =  x  +  y2. 


From  the  preceding  results  and  the  definition  of  r  Eq.  (2.8),  we  have 


ro=XUo  (~ 
\  v 


C  U0  \* 


Vx  / 


(3.12) 


We  can  use  Eq,  (3.12)  to  determine  the  skin  friction  when  U1  ^  0,  by  expanding 
rQ  in  a  Taylor  series  about  U0.  We  have,  if  T(U0  +  UJ  denotes  the  skin  friction 


+fhis  result  can  be  obtained  from  Ea.  (21)  of  [  1] It  also  was  given  explicitly 
by  Lagestrom,  Cole  and  Trilling  [  4J . 
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for  a  flow  U0  +  Ut  at  infinity  that 


r(U0  +  U,)=  T0  +  Ut 


This  result  is  correct  up  to  quadrate  terms  in  Uj ,  and  since  this  is  the  same  as  the 
accuracy  in  Eq.  (2.6)  we  may  use  this  result  to  compute  when  u>  =  0.  We  have 


=0) 


i  «,<*,» -0)-| 


( 


n\>x  J 


(3.13) 


These  results  will  be  of  interest  in  the  discussion  in  Sec.  5. 
4.  Case  Un  =  0 


In  this  section  we  require  that  Ut  «  (w  vJ)  A  as  indicated  just  after  Eqs.  (2.4). 
When  U0  =  0,  Eq.  (2.6)  reduces  to 


AA  f-  /32  A  ^  =  0  . 


(4.1) 


We  should  note  that  this  equation  gives  the  first  approximation  to  the  time 

j/ 

dependent  flow  even  when  U0  o  if  the  parameter  <r  =  (cov)  VcU0  is  large. 

This  can  be  seen  by  introducing  appropriate  dimensionless  variables  in  Eqs.  (2.5) 
and  (2.6)  and  then  expanding  ~f  in  powers  of  <r .  The  length  scale  that  one  should 
use  is  (V/w)1. 

To  solve  (4.1)  we  proceed  in  exactly  the  same  manner  as  described  in 
Sec.  3.  The  Fourier  transform  of  Eq.  (4.1)  is 

(|2  +  f  +/32)  (^2  +  n2  +k2)t(?  ,v  =  IV  f|0(?)  ■  (4.2) 


We  have  written  f10(f )  instead  of  fx  (T)  to  indicate  that  this  is  the  special  case 
UQ  =  0,  (i.e.  the  limit  case  <r->  co).  From  Eq.  (4.2)  we  obtain 

?  y  exp  W*) 

T' s>  J/  I  y  I  (£2  .  Jj2j  o  '  ' 
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and  finally 


f,o(f)  _  N+(1«)N.(?) 

2  =  i(f  -  ia)  * 


(4.4) 


where  N(f)  =  (f 2  +  /S2)‘i  +  (£2  +k*)%  is  analytic  in  the  strip  -  k*Im(?)  <-k. 

Before  we  proceed  we  must  determine  N.(1).  In  the  previous  section 
we  could  determine  K„(  f)  by  inspection;  here  the  situation  is  more  difficult. 
The  details  of  the  computation  of  N.I'?)  are  presented  in  Appendix  A,  It  is 
found  that 


r 


N.U)  =  exp 
where  in  the  limit  as  k  - 

F(?)  = 


2iri 


FO)  d? 


(4.5) 


1 


(Y+P)K 


£n 


/3  -  (**+*•>* 


(4.6) 


The  behaviour  of  F(?)  in  the  complex  plane  is  quite  complicated.  The  points 
£  =  i/3  and  ^  =  0  are  branch  points  of  the  function;  however  the  behaviour 
of  the  logarithmic  term  in  (4.6)  must  be  such  that  ^  =  -  i/3  is  not  a  branch 
point.  (Recall  N_("£)  must  be  analytic  in  the  LHP.)  Thus  we  choose  the 
branch  of  the  logarithm  by  defining  the  logarithmic  term  in  (4.6)  to  be  zero 
at  ^  =-i/3,  in  which  case  it  becomes  -  2ffi  at  ?  =  +  i/3.  This  can  be  seen  by 
writing  F(|)  =  2(^2  +  /32) jj/3  +  (?2  +/32)'^]/e7ri//2fj-  where  we  have 
taken  (-1)  =  e  .As  our  branch  cut  we  shall  take  the  entire  positive  imaginary 
axis.  The  determination  of  f  and  ("S  2  +  fi2)'*'  on  the  branch  cuts  is  given  in 
Appendix  B  where  the  results  are  used  in  evaluating  f10(x).  Finally  from  Eq.  (4.5) 
it  is  clear  that  N„(o)  =  1,  and  since  when  k  — >  0,  N(o)  =  j8  we  see  that  N+(o)  =  /3. 

We  may  now  compute  f10(x).  Upon  substituting  in  Eq.  (4.4)  and  using 
the  usual  inversion  formula  we  have,  letting  k  — ►  0,  a  — >  0  that 

i  tx  r  rf 


1  f,n(x)  = 

2  10  2  n 


exp 


1 

2ir  1 


F(?)d|  , 


(4.7) 
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where  the  path  of  integration  must  pass  below  the  origin.  This  integral  is 
evaluated  in  Appendix  B.  We  obtain 


l  l,0(x)  =  £  U  +  I  e’^  H(y)  dy  +  !  ,***'  Ply)  dy 


Here  x'  =  /3x  =  (1  +  i)  (o>/2V)Ax, 


(4.8) 


H(y)  =  exp 


1 

7T 


t  In  cot  -+  f — §— 
2  J  sin© 
0 


P(y)  =  exp  {  -  - 
n 


6 


sin  0 


L 


} 


t  =  arc  sin  y  , 
m  =  arc  sin  (1/y)  , 


and  0  <  y  <  1  implies  0<t<  jt/2;  also  1  *  y  <•  oo  implies  tt/2  >m>0. 

It  is  interesting  to  note,  using  Eq.  (4.8)  that  fiG(x)  is  of  the  form  0 
times  a  function  of  /3x.  The  expression  for  f10(x)  given  by  Eq.  (4.8)  may  be 
used  for  numerical  computation  when  (a>/2V)*-x  is  neither  large  nor  small. 

Of  course  the  second  integral  must  be  converted  into  an  integration  over  a 
finite  range  but  this  is  not  a  difficult  matter.  Using  Eq.  (4.8),  f10(x)/2  has 
been  computed  for  hx  =  1/2  and  1,  where  h  =  (w/2V)'i  .  The  results  of  this 
computation  are  given  at  the  end  of  this  section.  When  x  =  0  it  is  clear  that 
the  second  integral  in  (4.8)  does  not  converge;  this  of  course  is  not  surprising 
since  we  could  anticipate  from  the  result  of  Sec.  3  that  fl0(x)  should  behave 

's 

like  x’*  near  the  leading  edge  of  the  plate. 

To  complete  our  investigation  of  this  case  we  shall  obtain  expressions 
for  the  behaviour  of  f1Q(x)  for  large  and  small  values  of  hx .  To  do  this  we 
first  note,  since  f10(x)  is  identically  zero  for  x  0,  that  the  Fourier  transform 
of  flQ(x)  is  essentially  the  usual  Laplace  transform  but  rotated  through  90°. 
Hence  we  may  make  use  of  the  Tauberian  theorems  relating  the  behaviour  of 
the  Laplace  inverse  at  the  origin  to  the  behaviour  of  the  function  at  infinity, 

and  the  behaviour  of  the  Laplace  inverse  at  infinity  to  the  behaviour  of  the 
function  at  the  origin.  (See  [  5] .)  Thus  we  must  study  the  behaviour  of  N.(| ) 
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for  large  ^  and  small  £  . 

The  expression  for  N„(  £)  given  by  Eq.  (4.5)  can  be  written  in  the 
more  convenient  form  w 

N.(1)=  js  +  (1  +  s2)*!  exp  j-  j/n(coth  0/2)  dS^j  >  (4.9) 

where  s  =  ?//3  and  6  =  sintf's.  This  result  is  obtained  by  rationalizing 
{/3  +  (| 2  +j82)Ji}/|)3  -  (£2  +/32)*}  in  (4.6),  and  then  making  the  substitution 
s  =  sinh  6  .  it  is  not  difficult  to  show  from  (4.9)  that  as  s  — >  0 

N.(£) — *  1  j  (4.10a) 


and  that  as  s 


(^l*)  e’i/4  +  ot<^/5  ■  (4,10b) 


Using  these  results  and  Eq.  (4.4)  we  obtain  that  as  x  — ►  oo 

|  f  10(x)  —*■  P  -  (zf  enl/i  , 


and  as  x- 


In  particular  for  very  small  /3x 


\  'lo(X)  =  (M/  f1+  IT  +0^x)2] 


(4.11a) 


(4.11b) 


I  ‘.oM  -  \{ff  i'- 


*  e7Ti/B 


As  could  have  been  anticipated  Eq.  (4.11a)  is  the  result  for  an  infinite 
lamina  undergoing  simple  harmonic  motion.  (See  [  6] .)  This  result  represents 
a  wave  of  transversal  vibration,  i.e.  a  shear  wave.  Exactly  how  far  down  the 
plate  it  is  necessary  to  go  before  (4.11a)  is  valid  can  be  obtained  by  evaluating 
fl0(x)  from  Eq.  (4.8).  The  argument  and  a  dimensionless  modulus  of  f10(x)  have 

been  computed  using  Eqs.  (4.8)  and  (4.11b).  The  results  are  tabulated  in  the 
table  at  the  end  of  this  section,  and  are  exhibited  graphically  in  Figs.  3  and  4. 

It  appears  from  Fig.  3  that  Eq.  (4.11a)  is  an  excellent  representation  for  f10(x) 
for  hx  >  2.5  approximately.  Also  from  the  tabulated  results  it  appears  that 
Eq.  (4.11b)  is  valid  for  hx  *  1/2  approximately. 
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To  summarize  our  results  we  have,  under  the  assumptions 
(u>  V)*  »U0  and  (wvf>?  Ut,  that  the  skin  friction  t  is 

«•  -  /  U0  ^  +  /m,  {  fl0(x)  +  0  (•!,)}  e^4  (4.12) 

where  f0(x)/2  is  given  by  Eq.  (3.11)  and  flQ(x)/2  is  given  by  Eq.  (4.8).  Near 
the  leading  edge  of  the  plate,  t  takes  the  form 


r  =  f  V*.  {^)i(f-)fel(“t  +  V8)  •  <4-13> 


Far  down  the  plate,  i.e.  hx  >  2.5  approximately, 

r=  AUo  ♦/».  (ff  e1(“t  + 


(4.14) 


Thus  the  flow  consists  of  the  sum  of  a  mean  flow  plus  an  oscillatory  motion 
which  is  independent  of  the  mean  motion.  The  phase  advance  in  the  skin 
friction  increases  monotonically  from  n/ 8  at  the  leading  edge  to  jt/4  in  a 


distance  of  hx  =  2.5  approximately. 


TABLE 


Arg  ff l0(x)}  |  f l0(x)|  /2  (w/2  V)  * 


Evaluated  by 


Asymptotic  Formula  (Eq.  (4.11b) 
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5.  Case  Un  +  U,  e1^* 

Now  we  are  concerned  with  solving  Eq.  (2.6)  without  any  restriction 
on  the  parameters  involved.  In  theory  we  could,  of  course,  proceed  exactly 
as  we  have  in  the  past  two  sections.  However  because  of  the  non-homogeneous 
term  as  well  as  the  more  complicated  differential  operator  which  is  applied 
to  "f  we  should  expect  the  details  of  the  computation  to  be  more  difficult.  In 
actual  fact  a  direct  application  of  the  technique  used  in  the  previous  sections 
does  present  some  rather  formidable  problems;  hence  it  appears  desirable 
to  attack  this  problem  in  a  slightly  different  manner. 

First  we  note  that  'f'  can  be  expressed  as  the  sum  of  an  appropriate 
solution  of  the  homogeneous  equation  plus  a  solution  of  the  non-homogeneous 
equation.  Hence  we  can  write 

t,(*,y)+ r,{*,y)  ,  (5.1) 

where  ^  is  a  solution  of  the  homogeneous  equation  which  also  satisfies 
the  boundary  condition  ^^/dy  =  -  1  on  the  plate.  Then  the  solution, 
of  the  non-homogeneous  equation  (2.6),  must  satisfy  the  condition  y  -  0 

on  the  plate.  The  discontinuity  condition  in  'f'  can  be  expressed  as 

yy 

f,(x)  =  flt(x)  +f12(x)  where 

^yy(x,  o+)  =  -  V;yy(x,  o  -)  =  |-  fu(x)  , 

% yy(x,  0+)  =  -  )^2yy(x,  0  -)  =  ^  f12(x)  . 

First  we  shall  treat  ^  •  In  determining  fu  (x)  we  shall  proceed 
exactly  as  before.  The  Fourier  transform  of  is 

<12  +  !2  +  lt2){n2+?2+laW2}  tu ,  1)  ‘  .  (5.2) 
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Solving  for  Y^Jl)  and  inverting  over  we  obtain 

y.v  .  y  exp{-lyl(;2-t-k2)}  -exp{-|y[({2+  ia?+£2)*-}  Fu (| ) 
* ‘5,y)='iyl  (ia?+^-k2)  2 


7 


(5,3) 


Then 

*n(?)  _  M+(ia)M.(f) 
2  i(  |  -  ia) 


(5.4) 


where  M(£ )  =  (f 2  +  iaf  +  02)“*  +  (tf2  +  k2)*  .  The  splitting  of  M(|)  can  be 
carried  out  in  the  same  manner  as  was  done  in  Sec.  4,  hence 


JLn  ■  (z2  +  iaz  +02)^  +  (z2  +  k2)^  • 

- k - - - L  dz  . 

z  -i 

i 

If  we  multiply  the  top  and  bottom  of  Eq.  (5.4)  by  N+(o)  N_(  £),  note  that 
M+(o)  =  M(o)/M„(o)  =  0M.(o),  N+(o)  =  0,  and  N.(o)  =  1,  and  use  Eq.  (4.4)  we 
may  express  f u  ( J )  as 


/nM-<5>=  aJT 


(5.5) 


full) 

2 


f  M.(l)/N.(?)  f,n($) 

|_M.(o)/N.(o)  j  2 


(5-6) 


Thus  the  transform  of  the  skin  friction  contribution  from  the  homogeneous 
solution  can  be  expressed  as  a  function  of  t  times  the  transform  of  the  skin 
friction  when  UQ  =  0.  Equation  (5.6)  will  be  examined  in  more  detail  shortly, 
but  now  let  us  turn  our  attention  to  determining  the  non -homogeneous  solution 
of  Eq.  (2.6). 

Using  Eq.  (2.5)  it  is  clear  that  -(a/02)  is  a  solution  of  Eq.  (2.6). 

However  an  inspection  of  the  properties  of  X  as  given  in  Sec.  3  shows  that 
for  y  =  0,  Yz  =  -  (a/02)  X  will  lead  to  a  velocity  which  becomes  infinite  at 
the  leading  edge  of  the  plate,  i.e.  Xw(x>  o)  behaves  like  x‘*  for  small  x  >  0. 

This,  of  course,  we  can  not  accept.  But  we  can  add  solutions  of  the  homogeneous 
equation  to  this  solution,  and  the  result  will  still  be  a  solution  of  the  non-homogeneous 
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equation.  Hence  we  might  be  able  to  construct  an  appropriate  non -homogeneous 
solution  by  finding  a  homogeneous  solution  which  will  destroy  the  undesirable 
singularity  in  X  .  Thus  we  take  ^  to  be  of  the  form 


^2(x,  y)  = 


(x,  y)  +  </>  (x,  y) } 


(5.7) 


where  <f>  (x,  y)  must  satisfy  the  homogeneous  equation  (A  -  a.£/Zx  -  /32)  A0  =  0 
and  also  remove  the  singularity  in  X  .  To  describe  our  problem  completely 
let  us  specify  exactly  the  conditions  which  we  must  require  of  )^2.  The  first 
condition  is  V^2(x>  °)  =  0.  Secondly  we  require  that  y  =  0  on  y  =  0, 

x  >  0,  since  ^ir1/d  y  =  -  1  there.  Finally  <)2  y2/<^  y 2  =  Oony  =  0,  x  /.0, 
since  ft  (x)  =  0  for  x  <•  0.  And  of  course  also  '/"2 1  s  behaviour  is  such  that  its 
contribution  to  the  velocity  is  finite  at  all  points.  The  first  three  of  these 
requirements  are  clearly  satisfied  by  Xx  but  not  the  last  as  mentioned  earlier. 
Hence  <f>  must  satisfy  the  first  three  requirements  and  also  fix  the  singularity 
in  Xx.  The  determination  of  (x,  y)  and  in  particular  <L  (§,  o)  is  given  in 
Appendix  C . 

Now  the  transform  of  fl2(x)/2  is  given  by 

a  f  —  1 

~T“  =  ”F  ,y  =  0)+*yy($,y  =  0)J  ,  (5.8) 

and  noting  that  X^yd  ,  y  =  0)  =  i  f  f"0(  f  )/2,and  substituting  for  £yy(  ?  ,  y  »  0) 
from  Appendix  C  we  have 


fiaU.)  _  a 

2  '  “  Y 

where  a,  =  |(a2  +  4/32)^ 
and 


■  1 1  hm . 
a  2 

-  a}/2,  Q(5  )  =  {(J2 


a*  (S-ia,)*  Q-(  I  )|  ,  (5.9) 

+  ia  5  +J92)1  +  (52  +  k2)*)/{2(?2  +  ia 


Q.(?)  =  exp  (-  j}r  f  =1^  dx  ■ 
^  -06  ^ 


,  Im(^)  -t  0  . 
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Combining  Eqs.  (5.9)  and  (5.7)  we  see  that 

£,(?)  fM.(t)/N.a) 


M.(o)/N. 


(5.10) 


We  might  recall  in  considering  this  equation  that  fQ(\  )/2  is  the  skin  friction  when 
Uj  =  0,  and  f10(£  )/2  is  the  skin  friction  when  UQ  =  0,  (wV)^  >  Ut. 

To  invert  Eq.  (5.10)  in  general  is  a  rather  difficult  task;  however  it  is 
possible  to  obtain  expressions  for  ft(x)  for  small  x  and  for  large  x.  First  let 
us  consider  fu(£  )/2.  It  is  convenient  to  use  the  expression  for  fu  (?)/2  given 
in  Eq.  (5.6).  As  £  — 0,  (i.e.  x  — ►  oo)  the  bracketed  quantity  in  Eq.  (5.6)  is 
invariant  with  changes  in  a,  and  as  a  matter  of  fact  approaches  the  value  one. 

Here  the  first  approximation  to  fu(x)/2  for  large  x  is  simply 


t  n  (x) 
2 


A/ 


fio(x) 


ffi/4 


as  x 


oo. 


On  the  other  hand  M  .( I  )/N  .( $ ) 


1  as  $ 


(5.11) 

oo  and  hence  we  have 


*n  (x)  ^  N.(o) 
2  ~M.(o) 


N(o) 

/v— 

M.(o) 


*to(x) 


,  if 

(A  ( 

-)  I 

7TX  \ 

V 

y)  J 

ni/8 


as  x 


0  . 


(5.12) 


The  ratio  N„(o)/M_(o)  is  a  function  of  the  parameter  <r2  =  wV/ (c  U0)2.  It 
is  computed  in  Appendix  D  where  we  obtain  upon  letting  N.(o)/M.(o)  =  exp^I(<r)} 
that  for  <r  £  1/2 

i(,j>)  =  -  j_  (ifeL 

K  ]  it  Zl  (2n  + 1)! 


and  for  o' *  1/2 


2n+  1 
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where  X  =  2’1  e77^.  From  these  expressions  for  1(0  we  can  already  see 
the  effect  of  the  mean  flow  on  the  oscillatory  motion.  For  instance  as  U0  —  >  co, 

(i.e.  <y — >  0),  1(0  - >  *[■*  Jri/8  -  2‘  i,n  (20}  and  uP°n  substituting  in 

Eq.  (5.12)  we  see  that  the  phase  advanced  in  fu(x)/2  is  cancelled.  We  shall 
discuss  this  in  more  detail  after  we  have  obtained  appropriate  representations 
for  f12(x).  The  real  and  imaginary  parts  of  1(0  are  given  graphically  in  Fig.  5. 

Now  f12(0  is  given  by  Eq.  (5.9).  To  determine  il2(x)  for  small  x  we 
must  know  the  behaviour  of  f12(  I)  as  E — *«oo.  In  particular  we  must  study 
Q.(l)  for  £  large.  This  is  done  in  Appendix  E  where  we  obtain  the  result  that 

«*-«>■  x*  It  i  +0(f)  (5-13) 


where 


=  -ia, 


tr 


jri  +  1  + 


f  2  l-&\  ,  1 

U-r  2ir*J log  i  -  u*  j 


(5.14) 


and  S  =  aj/a2,  a2  =  {(a2  +  4  |32)*  +  a}/2.  Upon  substituting  for  Q_C§)  and  fQ(  |) 
in  Eq.  (5.9)  we  obtain  that  for  large  £ 


k 

a^-  2il 

2pz 


(5.15) 


where  G'  is  the  bracketed  quantity  in  Eq.  (5.14).  Finally  using  the  definitions 
of  a.lt  a2,  $  ,  £2  and  (T  Eq.  (5.15)  can  be  written  as 


H(tr') 

2 


+  o«'*) , 


(5.16) 


where 


H(c*)  =  1  +  1|  21  - 

n  2  6 


log  HtS  +i<1  'S) 


and  £  =  1  +  i  {(1  +  4i<r2)Ji  -  l}/2o-2.  Consequently 

as  x  — 

2  2  V.  'fvx/ 


(5.17) 
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We  are  now  in  a  position  to  discuss  the  behaviour  of  fl(x)/2  for  small 
x.  We  have  from  Eqs.  (5.12)  and  (5.17)  that  for  x  small 

¥'",w  •¥($#«■*>. 

(5.i8) 


where 


D(<r)  =  +  (2<r)*  e1^  e^1^8 


Thus  near  the  leading  edge  of  the  flat  plate  the  total  skin  friction  is 

iojt 


T-S-VaflUo)*  +/^u,D(r)  fSHa)* 

V  ttVx  /  V  / 


(5.19) 


As  U0 — v  cd,  or  w - >0  the  parameter  o' — y  0  and  D(<r)-^{(3/2  -  7cr2/127r)+  i^/ir2} 

and  hence  when  <?  =  0  (Le.  u>  =  0) 


-^(t ’fj  (1+  !  it) 


(5.20) 


This  result  is  in  agreement  with  Eq.  (3.13).  As  a  matter  of  fact  it  will  be  shewn, 
at  the  end  of  this  section,  using  Eq.  (5.10)  that  the  above  result  is  valid  for  all 
x  when  <?=  0. 

On  the  other  hand  when  cr1  — y-  oo,  i.e.  u>  — >  oo  we  have 
D(<T)  — *  (2£T)*  {l  +  (1  -  i)/(2ir <r)}  e 171/8  and  hence 


r=/*U, 


cU, 


-f 


jtVx  / 


ru. 


i(u>  t  +  71/ 8) 


(5.21) 


This  is  in  agreement  with  Eq.  (4.12)  as  could  have  been  anticipated.  As 
increases  from  0  to  co,  arg{D(<0}  and  hence  the  phase  advance  in  the  fluctuating 
skin  at  the  leading  edge  increases  monotonically  from  0  to  jt/8.  This  result  is 
presented  graphically  in  Fig.  7.  Also  the  modulus  of  D(<r)  is  given  in  Fig.  8. 

It  can  be  seen  from  Fig.  8  that  the  modulus  decreases  for  a  short  range  of  <T 


t 

I 
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as  o'  increases.  This  decrease  ie  extremely  slight,  of  the  order 
of  two  or  three  in  one  hundred  and  fifty.  However  the  asymptotic  formula  for 
DfO  for  small  o'  bears  out  the  numerical  calculations.  There  is  apparently 
no  obvious  physical  reason  why  the  amplitude  of  the  skin  friction  should 
decrease  slightly  as  increases  from  0  to  roughly  1/2. 

It  remains  now  to  consider  the  large  x  behaviour  of  f12(x).  (We  have 
already  shown  that  fu(x)/2— ►  (o>/9 )  e7^4  as  x  — oo,  see  Eq.  (5.11).  So 
we  must  now  consider  the  small  I  behaviour  of  f12(|).  This  is  a  fairly 
laborious  task  and  is  perhaps  not  worth  a  detailed  analysis.  It  can  be  seen 
from  the  definition  of  Q.(£)  that  Q.(o)  is  finite  and  hence  (S  -  iaj*  Q_Os) 
behaves  like  for  small  £  at  worst.  Also  it  f0(£)/2  behaves  like  (i£)5 
so  f12(  £)/2,  (see  Eq.  5.9),  will  consist  of  terms  of  the  type  (i|)^  ,  (i  £)^  , 
and  so  on,  at  worst.  Hence  fu  (x)/2  will  be  made  up  of  terms  like  x  x  ^ 
for  x  large.  Consequently  as  long  as  a,  /  0,  (i.e.  u>  ^  0)  these  terms  should 
die  out  more  quickly  than  those  contributed  by  fu  (x);  hence 

ft(x) — *"  (yj  ,  as  x — >oo  (5.22) 

This  result  is  to  be  expected  on  physical  grounds  since,  as  long  as  r'Z  is 
bounded  away  from  zero,  far  enough  down  the  plate  the  flow  should  consist  of 
a  shear  wave  superimposed  on  the  mean  flow.  As  <r2  increases,  (i.e.  as  u> 
increases  or  UQ  decreases)  the  shorter  the  distance  down  the  plate  we  must 
go  in  order  for  (5.22)  to  be  valid.  In  the  limitmg  case  UQ  =  0  this  distance 
is  approximately  given  by  (u>/2  v3)^  x  =  2.5. 

Finally  we  shall  investigate  the  case  =  0.  As  has  been  pointed 
out  earlier  (see  Sec.  3)  we  expect  in  this  limiting  case  that  f,(x)  =  (3/2)  (cUq/ttvx)^. 
If  we  set  u>  =  0  and  hence  /32  =  0  in  M(  I)  (see  the  definition  of  M(£)  immediately 
following  Eq.  (5.4)  we  see  that  M(^)  becomes  identical  with  K(^)  and  hence  upon 
examining  Eqs.  (5.4)  and  (3.9)  and  (3.10)  we  obtain 
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(5.23) 


To  compute  f12(x)/2  when  02  =  0  we  turn  our  attention  to  Eq.  (5.9).  It  is 
simple  to  show  by  an  examination  of  Q(f )  and  R.(^  )  (defined  in  Appendix  C) 
that  when  02  =  0  and  hence  ax  =  0  that  the  bracketed  quantity  in  Eq.  (5.9) 
vanishes.  Hence 


lim 
— *0 


*!,(*) 

2 


*»(•*)  j 


This  computation  is  straightforward  and  we  obtain 

*»(*)  a* 

2  '  2(i$)* 

when  /32  =  0.  Consequently 

fi2(x)  1  /cU0\£ 

2  2  \  )  » 


(5.24) 


(5.25) 


and  adding  Eqs.  (5.23)  and  (5.24)  we  obtained  the  desired  result 

M*)  =  3  /  cU0 

2  2  ^  ttVx 

when  w  =  0. 

To  summarize  briefly  we  have  found  that  when  (wv*)^?*  UQ  and  U1  the 
mean  motion  and  the  oscillatory  motion  are  independent.  The  total  skin 
friction  is  given  by  Eqs.  (4.12)  and  (4.8).  Asymptotic  expressions  valid  near 
the  leading  edge  and  far  down  the  plate  are  given  by  Eqs.  (4.13)  and  (4.14). 

It  is  found  that  the  phase  advance  in  the  time  dependent  skin  friction  increases 
from  tt/8  at  the  leading  edge  to  jt/4  in  a  distance  of  approximately  {uj/2'\>)x  x=  2.5. 
Finally  when  the  interaction  of  the  mean  flow  and  the  time  dependent  flow  can 
not  be  neglected  a  result  valid  near  the  leading  edge  and  one  valid  far  down¬ 
stream  are  obtained  for  the  oscillating  skin  friction.  These  results  are  given 
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in  Eqs.  (5.19),  (5.20),  (5.21)  and  (5.22).  It  is  found  that  the  phase  advance 

in  the  time  dependent  skin  friction  at  the  leading  edge  increases  from  0  to  n/ 8 

2 

as  the  parameter  c  Uq/(o>  V)  increases  from  0  to  oo.  A  graph  of  the  phase 
advance  versus  this  parameter  is  given  in  Fig.  7.  Far  down  the  plate  one 
obtains  the  usual  shear  flow  solution  (see  Eq.  (5.22)). 


Appendix  A 
Computation  of  N .(f) 


In  the  strip  -  k  *Im  (§)  <  k,  N(?)  is  analytic  and  non -vanishing. 

Hence  using  Cauchy’ s  Integral  Theorem  we  may  write 

•^n  N(J)  =  ^h[n+(?).N.(?)7.  JL  fin  dz  ,  (A.l) 

n  c  z 

where  C  is  the  complete  contour  mdicated  in  Fig.  (1),  (i.e.  C=  f’,  +  Lx  +  C2  +  L2). 
If  we  allow  the  ends  of  this  contour  to  approach  infinity  at  an  equal  rate,  the 
integrals  over  and  f2  will  converge,  and  the  contributions  from  the 
integrations  over  Lt  and  L2  will  vanish  giving 

ZnN+(?)+/nN  (I)  =  ~  [***!(?}*£  +  JL  illZM  dz  .  (A. 2) 

t  2tti  i  •»”  **-s  r 


Jr 


z  -  l 


2ni 


z  -  f 


The  first  integral  in  the  above  equation  is  uniformly  convergent  for  Im(  ?)  >  -  1 
and  hence  represents  an  analytic  function  of  J  there.  A  similar  statement  is 
true  for  the  second  integral.  Associating  these  integrals  with  N+(^)  and  N^) 
respectively,  we  have 

r 


Jin  N.(5)  = 


2iri 


Xn  (z2  +  j32)^  +  (z2  +  k2)*-  J 
z  -  i 


dz 


(A.3) 


If  we  confine  ^  so  that  Im(f )  <•  0  we  can  take  P2  to  be  the  real  axis 
in  evaluating  (A.3).  Let 

F(?)  . 

2m  d  |  ’ 


then 


r  tg.±wi 

J.  (x-f)2 


<*-?)' 

Upon  integrating  by  parts  and  writing  this  as  the  sum  of  two  integrals  from  -  oo 
to  0,  and  0  to  oo  and  then  making  the  usual  change  of  variables  we  obtain 

2x2  dx 


F(?)  = 


r 


(A.4) 


A2 


In  the  limit  as  k  — >  0  the  value  of  this  integral  is 


F(?)  - 


1 

(?2+02)* 


An 


g  +  (f  a+g*)* 

13  - 


(A.  5) 


Now  we  can  specify  N.(? )  precisely  by 

rf 

F(?)d?|.  (A.  6) 

0  J 

The  lower  limit  in  the  above  integral  could  be  chosen  differently  but  this 
merely  multiplies  N„(?)  by  a  constant;  the  reciprocal  of  this  constant  would 
then  multiply  N+(t). 


N  Jf)  =  exp 


Appendix  B 
Evaluation  of  f10(x)/2 


Substituting  in  Eq.  (4.4)  for  N_(f ),  and  using  the  usual  inversion 
formula  we  have  letting  k  — >  0,  a  — >  0  that 

I  -  £  1 1 e2rr  exp  {•  iir JF{?)  dJ}d?  -  (B-1) 

-•e  L  Jo  J 

where  the  path  of  integration  must  pass  below  the  origin.  Making  the  change 
of  variables  s  =  £//3,  x'  =  /Sx  and  noting  that  for  x  >  0  we  can  deform  our  path 
of  integration  to  an  integration  over  S  (see  Fig.  2)  leads  to 


As  indicated  in  Fig.  2,  the  contour  S  is  broken  up  into  the  four  parts 
Sj,  S  2>S3,Sy,  plus  the  enclosures  about  the  branch  points  at  s  =  i  and  the 
origin.  It  is  not  difficult  to  show,  that  as  the  radii  of  these  enclosures  are 
allowed  to  approach  zero,  that  there  are  no  contributions  due  to  the  integrations 
around  the  branch  point  s  =  i,  and  that  the  integration  around  the  origin 
contributes  2 rr.  On  the  Si  we  have  the  following: 


B2 


Sa:  s  =  e"3,ri/2  y;  for  y  >  1,  (1  -  y2)*  =  e’77^2  (y2  -  1)* 

“£  C  °° 

G(s)  =  -  i{y  +  (y2  -  1)*J  exp  {+  I 


i _  (y2  1 

1)*  /n  (y2  -  l)*  +  i  dyJ  ’ 


S2:  s  =  e*3'1/2y,y<l 


G(s)  =  exp 


.  i.  4>[l  +  (l_-.ys)*]  dy  -  i  arc  sin  y  k 

.  1  J„  (1  -  y2)*  J 


_  ffi/2  . 

S3:  s  =  e  y,  y  <  1 


G(s)  =  exp 


1  ft n[l  +  (l  -y2)*]  -Ay 


(1  -y2)* 


dy  +  i  arc  sin  y  f  , 


S4 :  s  =  e 771/2  y,  for  y  >  1,  (1  *  y2)*  =  e 771/2  (y2  -  1)* 


G(s)  =  i-j  y  +  (y2  -l)1  V  exp 


’2  -  l)<  -  i 


2?ri  (y  -  1)^  (y2  -  l)*  +  i 


In  determining  these  expressions  for  G(s)  it  was  found  convenient  on  St  and  S¥ 
to  write  the  integral  0  to  y  as  the  sum  of  the  integrals  0  to  1  and  1  to  oo  minus 
the  integral  y  to  oo.  Combining  all  of  these  results  we  obtain  finally 

|  fio(x)=  \  +  J^e~yx'  H(y)  dy  +  lllJL l.lii^.1,5  P(y)dy|.  (B,4) 


Here 


H(y)  =  exp  4  -  - 


1  r  Jtn[l  +(1  -  u2)l  -inu 


(1  -  u2)* 


duL  y  <  1 


=  exp 


:-7[tincot  i  +/lI7e  d6]' 

.  Jo  * 


d© h-  ,  t  =  arc  sin  y  , 


and 


B3 


P(y)  =  exp 


(y*  •  D*  •  i 

(r  - 1)*  +  i 


,  y  >  i 


=  exp 


m  =  arc  sin  (1/y)  . 


When  y  ranges  from  0  to  1,  t  goes  from  0  to  jt/2;  similarly  l<y  ^cd  implies 
»r/2>m  5-0. 


Appendix  C 

Computation  of  y) 


The  x  Fourier  transform  of  a  function  0(x,  y)  which  satisfies  the 
differential  equation  (A  -  a  V<i  x  -  /32)  A0  =  0  and  the  condition  0(x,  o)  =  0  is 

•r,  y  expf-|y)(?2 +k2)’i} -exp{-|y|(S2  +  ia?  +  ,32)'4} 

* y  “  tf  ■  ,c-1) 


where  A(|)  is  an  arbitrary  function  of  I  .  The  condition  that  *%/Zy  =  0  on  y  =  0, 
x  >  0  requires  ({ ,  y  =  0)  be  analytic  in  the  UHP,  and  the  condition  d-2  f^y2  =  0  on 
y  =  0,  x  <  0  requires  0yy(  y  =  0)  be  analytic  in  the  LHP.  These  conditions  may  be 
written  as 


*y<?,  °»=  -  (g2+k2)„f  +  <1*-+ ij±t££  C  (^) 

y  A(S) 


(C.2) 


«yy(*,  °)  * 


(5  a -ka)  -(ta  +  lai  +  fi2)i  =D  (?) 
A(5) 


(C.3) 


where  C+(£)  and  D„(£)  are  functions  analytic  in  the  UHP  and  LHP  respectively. 
The  final  condition  on  may  be  expressed  as  X^x,  o)  +  0y(x,  o)  be  finite  at  the 
origin.  Since  the  behaviour  of  the  transform  at  infinity  determines  the  behaviour 
of  the  function  at  the  origin  (see  Sec.  4  and  [5]),  and  since  differentiation  with 
respect  to  x  implies  multiplication  of  the  transform  by  i  'S  we  may  express  this 
condition  as 


i|  Xy(x,  o)  +  <fy(x,  o) 


■(1+e) 


as  £ — *■  oo,  where  e  >  0.  Using  Eqs.  (3.5),  (3.10),  and  C.2)  this  becomes 

{e'v ft*  +c+<*>]  ~r(1*c)  -  <c-4> 

as  {—>oo.  Hence  C+(»~{’*  +  0[f'(1+e)]  as  $  — oo. 


C2 


We  can  determine  the  functions  C+('§)  and  D.(  I)  by  an  application 
of  the  Wiener -Hopf  technique.  First  dividing  (C.3)  by  (C.2)  gives 


D.(£)/C+(S)  =  -  {(S2  +  iaJ  +  02)*  +(|2+k2)i}, 


-  -  M(T)  , 


(C.5) 


where  M(?)  is  analytic  in  the  strip  -  k  <■  Im(t)  <k.  It  is  convenient  to  write 
(C.5)  as 


D.(1)/C+<!)  =  -  Q(*)R<?)  , 


(C.6) 


where  R(?)  =  2(£2  +  ia  ?  +  02)^  and  Q(f)  =  M(if  )/R($).  Now  we  can  split 
Q(£)  and  R(  f)  into  the  product  of  factors  analytic  in  the  UHP  and  LHP,  and  as 


a  result  we  obtain 


Q.(??)R-(|)  =  -C+(|)Q+(^)R+(^)  =  P(*) 


(C.7) 


where  P()j)  is  an  entire  function.  The  factors  R-C|),  R+(£),  and  Q-(£)  are 


given  by 


R-(?)  =  (?-ia,)*  , 

R+($)  *  2{g  +  ia,)1*  , 


(C.8) 


Q.(?)  *  exp  f-  2 JJ-  f  da  ■  for  Im(  V)  *0  . 

-o*  <- 


Here  at  =  £(a2+4/32)*  -  a}/2  and  a2  =  £(a2  +  4/32)^  +  aJ/2 .  To  determine 
P(^)  we  must  first  anticipate  that  Q_(£)  =  0(1)  as  t  — *-oo,  and  hence  since 
Q('5)  — +-  1  as  |  — ►  co  the  same  is  true  of  Q+Cf).  Also  since  D-(^)  is 
proportional  to  the  transform  of  the  skin  friction  we  expect  D.(^)-~<  0(^*0 
as  | — >oo.  Using  these  order  conditions  plus  Eqs.  (C.8)  and  the  information 
about  C+('|)  we  see  that  P($)  must  be  a  constant,  say  P1 .  Then 


C+(|)  = 


Qf(*)*M*> 


(C.9) 


,  D.($)-P,Q-(f)R..(*)  . 


C3 


To  determine  P1  we  substitute  the  above  expression  for  C+i'f)  in  Eq.  (C.4). 

-V  + 

Noting  that  R+(f )  ~  2|  and  Q.( })  =  1  plus  terms  of  higher  order  for 
large  £  it  is  clear  that  P'  =  a^  .  Consequently 

iyy(  ?  ,  o)  =  e'3”i/4  a*  Q.(?)R.(D  •  (C.10) 


As  was  pointed  out  earlier  we  anticipate  that  Q_(  ?)  is  0(1)  as  1  *co; 
and  since  Q-(  |)  is  determined  only  up  to  a  constant  multiplier  we  may 
take  the  leading  term  to  be  one.  Actually  the  expression  for  Q-U )  as 
given  in  (C.8)  does  lead  to  the  desired  expression  when  |  is  large. 
This  will  be  seen  in  Appendix  E. 


Appendix  D 
Evaluation  of  I(  rf 


The  functions  N .( |)  and  M.(|)  are  defined  by 

(in  f(z2  ♦|3‘)*  +  (z,  +  k,)*l 
|2?ri  1  z  -  ? 


N.(?  )  =  exp  < 


dz 


}• 


and 


M.($)  =  exp 


f  £n[(z2  +iaz  +/32)*  +  (z2+k2)* 

L  dz  ‘ 

\2ffi  * 

1  z  ’  1 

<■ 

If  we  restrict  £  so  that  Im(  t-)  <  0,  then 

ML>  «  exD  f+  J_  f.n  (x2  +  iax  +  g2)*+(x2  +  k2)*  te  1 
(  2iri  j  (x2  +  )92)'lt  +  (x2  +  k2)ii  x-? 

^  —at  — 

(D.l) 

In  the  limit  as  ^  — *  0,  our  integration  path  must  be  indented  to  go  above  the 
origin;  i.e.  the  path  of  integration  is  from  -  oo  to  -  e ,  then  a  semi -circle  of 
radius  c  in  the  UHP,  and  then  the  real  axis  from  e  to  oo.  In  the  limit  as 
e  — ►  o  the  integration  over  the  semi -circle  is  zero.  Letting  r  =  a x///3/2, 
r2  =  |/3 12/ a2  =  wV/(c  U0)2,  N„(o)/M_(o)  =  exp  I(<r),  and  noting  that  as  k  — *-0, 
(x2  +  k2)^  =  JxJ  we  obtain  from  (D.l) 


m  -  ;rr  f-tn  c  Llf..^r- 

J  [i(l  -  r)  +  <r2  r2]*  +<jr 


(D.2) 


This  integral  can  be  evaluated  by  first  computing  dl/dr  and  then  integrating 
with  respect  to  <r\  We  have 


OO 


[r2  r  2 


(D.3) 


( 


D2 


where  %  =  eir^*/2.  Consequently 


ee 


l(<r)  =  -i- 

2ii 


r 


xa  i*4 .  sr 

r  -  o  r 


(D.4) 


The  constant  of  integration  (i.e.  the  upper  limit  in  the  above  integral)  is 
determined  by  noting  that  when  a  =  0,  N  .(o)  =  M_(o)  hence  I(<r)  =  0,  but  a  =  0 
corresponds  to  rf  =  co. 

For  s'  1/2, 1(<r)  can  be  evaluated  directly  by  expanding  the  integrand 
in  a  power  series  in  X/t  and  then  integrating  termwise.  We  obtain 


®D 


KO  -  -  i  > 

*  £•  (2n  + 


1) 


2  > 


1/2  . 


(D.5) 


For  <T <.  1/2  the  evaluation  of  (D.4)  is  slightly  more  difficult.  It  is  convenient 
to  write  (D.4)  as  an  integral  from  IT  to  go  plus  an  integral  (r'  to  X ,  then  we 
obtain 


K<r)  =  -  \  in( 2r)+  te  [lintel  dr  . 

V 


(D.6) 


After  writing  the  integral  in  (D.6)  as  one  from  0  to  minus  an  integral  from 
0  to  o'  we  obtain 


i<o-  -  f ^ 


l  in  dr 

r  *  X  -  r 


v2n+  1 


tL  .  I  /n(2<r)  +i  _ 

8  2  .4  (2n  + 1)2  ’ 


<r«4/2  . 
(D.7) 


Appendix  E 

Asymptotic  Expansion  of  Q.(£) 


The  function  Q-(  |)  as  defined  by  (C.8)  is 


Q.(f)  =  exp 


1 

2ni 


in  Q(x) 

x-* 


dx  f  for  Im(S)  <  0  ,  (E.l) 


where 


Q(|)  = 


+  +  (fca+k*)* 

2(# 2  +  ia  *  + 


Except  for  branch  points  at  £  =  ik  and  ^  =  iax,  Q(|)  has  no  other  singularities 
in  the  UHP.  Hence  using  Cauchy' s  Integral  Theorem  the  integral  in  (E.l)  can  be 
converted  to  an  integral  over  T  (see  Fig.  6).  So  if  we  write  Q_(  \)  =  exp  {-G($)/2jri} 
we  have 


O(S).  -  dz  . 

For  large  £  ,  G(f)  has  asymptotic  representation  of  the  form 


(E.2) 


G(^  )  =  G0  +  Gx/£  + 


(E.3) 


where 


G0  =  lim  G(|) 

t 


Gi  = 


lim  V  ^ 

d£ 


(E  .4) 


It  is  not  difficult  to  show  from  (E.2)  that  G0  is  zero.  Performing  the  limiting 
operation  indicated  in  (E.4)  we  see  that 


r 

It  is  convenient  to  write  Q(z)  as 


V(z)  = 


(z  +  i  k)*  (z  -  i  k)^ 
(z  -  iax)*  (z+  ia^* 


(E.5) 


E2 


To  evaluate  (E.5)  in  the  limit  as  k — >-0,  we  make  the  change  of  variables 
x  =  z/iax ,  and  let  $  =  ax/a2.  Then  our  branch  cut  lies  on  the  real  axis 
from  0  to  1,  and  our  path  of  integration  is  clockwise.  The  function  V(z)  has 
the  values 


V  (x)  = 


_ x  e~^i/2 

(1-x)  (x  +  l/J)* 


V(x)  = 


X  e'”,i/2 

(l-x)*(jc  +  l/j)* 


> 


on  the  upper  and  lower  sides  of  the  branch  cut.  We  obtain 

i 

G,  =  -  ia,  f  in  tl-laW  dx;  q(x)  = 

J  1  -  iq(x) 


(l-x)*(x  +  l/f)’ 


(E.6) 


The  integral  in  Eq.  (E.6)  can  be  evaluated  by  integration  by  parts  and  some 
straightforward  algebraic  manipulation.  We  obtain 


{ 

1  *  & 


iri  +  1  + 


-in 


1  +  i  /  ' 
1-  i  i  * 


So  for  £  large 

Q-($) 


where  Gt  is  defined  above. 
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